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ABSTRACT 

Simultaneous measurements of fast electric field changes and RF radiation from first 
return strokes indicate that a time delay exists between the first noticeable RF radiation 
and the beginning of the return stroke as determined from the fast field change. Measure- 
ments in Florida showed this delay to be on the order of 20 ns in the HF-UHF range. 
Similar results have been found by this author in experiments in New Mexico (TRIP-79) 
and in Virginia at NASA's Wallops Flight Facility. 

It will be shown in this paper that the combination of the finite conductivity of the 
earth plus the propagation of the return stroke current up the channel results in an appa- 
rent time delay between the fast field changes and RF radiation for distant observers. 

The time delay predicted from model return strokes is on the order of 20 /zs and the 
received signal has the characteristics of the data observed in Virginia and Florida. 

The theoretical work has been done using a piece wise-linear model for the return 
stroke channel and a “transmission line” model for current propagation on each segment. 
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Radiation from each segment is calculated over a flat earth with finite conductivity using 
asymptotic approximations for the Sommerfeld integrals. The radiation at the observer 
is then “processed” by a model AM radio receiver. The output voltage has been calculated 
for several frequencies between HF-UHF assuming a system bandwidth (300 kHz) charac- 
teristic of the system used to collect data in Florida and Virginia. Comparison with the 
theoretical fast Held changes indicates a time delay of the order of 20 ps. 
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THE INFLUENCE OF GROUND CONDUCTIVITY ON THE STRUCTURE 
OF RF RADIATION FROM RETURN STROKES 


I. INTRODUCTION 

When measurements are made of RF radiation from the first return strokes of cloud-to- 
ground lightning flashes it is observed that the RF radiation does not start with the beginning of the 
return stroke as indicated by its fast electric field change but rather begins several 10’s of micro- 
seconds later. This time delay was first reported by Brook and Kitagawa (1964) and Takagi (1969a. 
b) in measurements made in New Mexico, and later was investigated by Le Vine and Krider (1977) 
in measurements made during the Thunderstorm International Research Project (TRIP) in Florida. 

Le Vine and Krider (1977) reported a relatively constant delay of about 20 ps at several frequencies 
between 3 MHz and 300 MHz. The existence of such a time delay is particularly surprising in return 
strokes because the electric field changes from first return strokes are known to begin quite abrupt- 
ly (Weidman and Krider, 1978;Lin et al., 1 979; Tiller, et al., 1976;Guo and Krider, 1982);conse- 
quently one would expect the beginning of the return stroke to be associated with radiation over a 
large range of frequencies, including RF. It is the purpose of this paper to show that this intuition 
is correct but that when the effects of the dispersive conductivity of the earth on the radiation from 
the return strokes is taken into account the effect is an apparent time delay. 

In the sections to follow, a review will be given first of the measurements establishing the exis- 
tence of the time delay. Examples will be presented from data collected by Le Vine and Krider in 
Florida and by Le Vine at Langmuir Laboratory in New Mexico and at NASA’s Wallops Flight Facility 
in Virginia. Then theory will be presented for radiation from return strokes in which the effects of 
a finitely conducting earth are included. The theoretical work employs a piecewise linear model 
for the return stroke channel and a “transmission line” model for the current on each segment. 
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This model is used to calculate the electric field radiated from the return stroke over a conducting 
earth using an asymptotic approximation which is valid for observers far from the channel. The 
electric field is assumed to be the input to an ideal AM radio receiver similar to those used in the 
experiments and a mathematical model for the receiver is used to calculate the receiver output vol- 
tage at several frequencies. Comparison of these calculations with the beginning of the return 
stroke indicates a time delay consistent with experiment. 
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II. THE TIME DELAY 


Figure 1 is a block diagram of an experiment of the type used by Le Vine and Krider (1977) 
and also by Le Vine (this paper) to observe the time delay. In these experiments measurements 
were made of the radiation from first return strokes of lightning flashes using a fast electric field 
change system (Uman, 1969; Krider et al., 1977) and also an AM radio receiver. The two signals 
v jre recorded simultaneously and stored on a digital sample and hold device (e.g., a Biomation 
model 8 1 00 waveform recorder) to obtain high time resolution. The analogue output of the wave- 
form recorder was displayed on a dual beam oscilloscope and photographed for a permanent record. 
Experiments of this sort have been performed by Le Vine and Krider in Florida (Le Vine and Krider. 
1977) and by Le Vine at Langmuir Laboratory in New Mexico during TRIP-79 and at NASA’s 
Wallops Flight Facility in Virginia during the summers of 1979 and 1980 as part of NASA’s Storm 
Hazards program (Pitts, 1982). The instrumentation used in all of these experiments was similar to 
that used by Le Vine and Krider (1977) in Florida. The RF receivers used at frequencies greater 
than 30 MHz were commercially available AM receivers (Watkins-Johnson models WJ-977 and WJ- 
S^O with WJ-9082 tuning heads) and at frequencies of 3 MHz and 30 MHz they were fixed tuned 
receivers specially designed by the Georgia Institute of Technology (Le Vine et al., 1976). The out- 
put of each receiver was filtered to provide a signal with a common bandwidth (300 kHz) at all fre- 
quencies. In the measurements in Florida, dipoles were used at all frequencies, but in the later 
measurements the dipoles were replaced with disk cone antennas at frequencies greater than 50 
MHz. The fast field change system was designed by Krider (Krider, et al., 1977) and consisted of a 
flat plate antenna about 30 cm in diameter followed by an integ/ator. The integrator compensates 
for the capacitive impedance of the electrically small antenna so that the output voltage is propor- 
tional to the incident electric field, E(t). The time constant of the integrator was about 1 ms, which 
determined the low frequency response of the system (about 1 kHz) and the high frequency re- 
sponse (several MHz) was determined by the operational amplifier used in the integrator. It is in this 
frequency range that most of the energy radiated from return strokes is to be found (Oh, 1969: 
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LeVine, 1980;Kimpara. 1965;Homer, 1 964; Cianos, Oetzel and Pierce, 1972). Consequently, the 
fast field change system provides a measure of the overall shape of the waveform radiated from re- 
turn strokes. The high frequency detail it lacks is the information provided by the radio receivers. 

Figures 2 and 3 are examples of lightning data collected in Florida with this system. The lower 
trace in each figure is the fast field change from a first return stroke and the upper trace is RF radia- 
tion at 295 MHz. These examples are part of the data set collected by Le Vine and Krider during 
TRIP-76 (Le Vine and Krider, 1977) but not previously published. The fast field change consists of 
a sequence of small pulses preceding an abrupt jump in the electric field (center). The large field 
change in the center of the record is characteristic of fast electric field changes associated with first 
return strokes (Weidman and Krider, 1978; Tiller et al., 1976; Lin et al., 1979) and the preceding 
sequence of pulses are typical of radiation from the stepped leaders which proceed the return stroke 
(Krider et al., 1977). The top trace in each figure is the voltage output of a radio receiver tuned to 
295 MHz. (Figure 3 is an example of horizontally polarized radiation. All other examples in this 
paper are for vertically polarized radiation.) Clearly, there is detectible radiation at this frequency 
during the return stroke. However, the radiation does not begin with the abrupt transition in the 
incident E(t) as one would expect; rather, it begins 10 or 20 ps after the initiation of the fast field 
change and does not reach a maximum (the horizontal lines indicate receiver saturation) for several 
10’s of ^s after the beginning of the return stroke. A similar delay between the beginning of the 
fast electric field change and the beginning of the RF radiation was observed at all frequencies (3- 
300 MHz) at which measurements were made in Florida (Le Vine and Krider, 1977). Similar results 
have also been obtained by Le Vine in measurements at NASA’s Wallops Flight Facility in Virginia 
and at Langmuir Laboratory in New Mexico. Figures 4 and 5 are examples of vertically polarized 
radiation collected at 3 MHz (Figure 4) a id at 139 MHz (Figure 5) in Virginia and Figures 6 and 7 
are examples of the fast elertric field change and RF radiation at 3 MHz and 170 MHz, respectively, 
collected during TRIP-79 at Langmuir Laboratory in New Mexico. The time delays observed in 
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these experiments were on the order of 20 ms except in Virginia where the delay at 1 39 MHz was 
closer to 40 ms. A consistent feature of the data is the rather abrupt beginning of the RF radiation 
at the higher frequencies (e.g., Figures 2, 3, and 5) but the gradual increase in signai level at 3 MHz 
(e.g., Figure 4;Le Vine and Krider, 1977, Figure 2). 
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III. RECEIVER RESPONSE TO RADIATION FROM RETURN STROKES 


It is the objective of the remainder of this paper to show that a time delay in the RF radiation 
from return strokes is a natural consequence of the dispersive dielectric constant of the earth and 
the propagation of the return stroke current up the lightning channel. This will be done theoreti- 
cally by adopting a model for the return stroke, using this moo '1 to compute the radiated electric 
field, and then processing this radiation through a model RF receiver similar to those used in the 
experiments. 


The standard AM radio receivers used during the experiments are devices which amplify and 
detect the envelope of an amplitude modulated sinusoid. In practice this is normally done by trans- 
lating the input signal to an intermediate frequency where the actual detection is done. However, 
the net effect is simply that of an envelope detector in series with a filter which represents the effec- 
tive bandwidth and gain of the system. It can be shown (Appendix A) that the voltage, V Q (t), out 
of such a system tuned to a frequency, v 0 , is 


V 0 (t) = 



a iv Q + £)E 2 (* 0 + f)H(f)e" j2,rft d? 


(1) 


where the bars 1 | indicate the magnitude of the complex number inside, E z (v) is the spectrum 
(Fourier transform) of the vertical component of the electric field radiated to the receiver from the 
return stroke, a(v) represents the effect of the antenna and system gain, and where H(i>) is a filter 
representing the equivalent bandwidth of the system. Equation 1 is obtained under the assumption 
that the bandwidth of the system is small compared to the nominal frequency, v Q , to which the re- 
ceiver is tuned. 


In order to compute the receiver output voltage, V 0 (t), the spectrum E z (tO of the incident 
radiation is required. This will be computed here by adopting a piecewise-linear version of what is 
commonly called a “transmission line” model for return strokes. The first step in this procedure is 
to create a channel out of linear segments placed end-to-end in a computer generated random walk. 
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The procedure is suggestive of the manner in which channels are formed in real lightning by the 
stepped leader. In the implementation of this procedure, a probability distribution is chosen for the 
change in cartesian coordinates that occur in going from the beginning of one segment to its other 
end. The coordinates themselves are selected by the computer “at random” from this distribution. 
Details and examples are given by Le Vine and Meneghini (1978a). A change in this procedure 
adopted by Le Vine since 1978 has been to use two probability distributions. The second distribu- 
tion is used to vary the statistics (the mean) of the first distribution. This gives a channel with small 
scale fluctuations (tortuosity) determined by the primary distribution and these are superimposed 
on large scale randomly changing drifts determined by the second distribution. Figure 8 is an exam- 
ple of a channel generated in this manner. The second distribution determines the rather long 
(several hundred meter) drift of the channel to the left or right. The two scale model was adopted 
after statistical comparison of computer generated channels with photographs of lightning channels. 
The two scale model not only yielded the most realistic looking channels, but it was also necessary 
to get the correlation between junction points to agree with the data (photographs). A two scale 
channel model was also found to be necessary by Ribner and Roy (1981) to achieve satisfactory 
simulation of thunder. 

Given the channel, the return stroke is modelled by a current pulse which propagates along 

each segment with constant velocity. That is, it is assumed that the current on each segment has the 
_ /\ 

form I(t - fij • rj /vp where fij is a unit vector along the i-th segment in the direction of current flow 
and Vj is the velocity of propagation of the current pulse along the i-th segment. In the literature on 
lightning, this model for the current is called a “transmission line” model. It was originally sugges- 
ted by Dennis and Pierce (1964) and developed formally by Uman and colleagues (Uman et al., 
1973; Uman and McLain, 1969; McLain and Uman, 1971) who also demonstrated agreement with 
data (Uman et al., 1973). The var cion employed here is a piecewise-linear adaption of the model 
which permits the effects of channel tortuosity to be included (Le Vine and Meneghini. 1978a; Le 
Vine, 1980; Le Vine and Meneghini. 1983). 


The fundamental analytic problem in this model is to calculate the electromagnetic fields radi- 
ated from each segment. Given the solution the computer can be programmed to obtain the total 
field at a particular observation poim by summing the contribution from each segment, keeping 
track of the phase as the current pulse propagates up the channel. This will be done here for an 
observation point near the ground and far from the channel compared to the channel length. A 
solution will be obtained for each element which takes into account the conductivity of the earth. 

The effect of the conductivity of the earth is included in the calculations by using Norton’s 
form for the solution for radiation from a dipole above a conducting earth (Norton, 1937). In the 
case of an observer very far from the element compared to the height of the element above the sur- 
face asymptotic approximations result in a particularly simple form (Appendix B) and the solution 
for the vertical (Le., z-component) of electric field at the surface can be written: 

E Z M => A h E“ z « + VQW (2) 

where Ej^(i>) and E^(iO are the solutions for the z-components of electric field due to a filament 
abo ve a perfectly conducting (o = <») ground plane oriented parallel (E^) or perpendicular (E^) to 
the surface. The A^ v are correction factors which account for the finite conductivity of the 
ground: 

Av - < 3a > 

a h=7 + 1 [ 1+j l^W] (3b) 

where f = (z'/RH/e^ and e r is the (complex) relative dielectric constant of the earth (e r = ? r ' + j e r "), 
z is the height of the element above the surface, R is the distance from the element to the observer 
and k = 2rri>/c. 
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ORIGINAL PACK r ,3 
OF POOR QUALiVY 


1 he E“ z (i>) and E*(y) in Equation 2 are the solutions for radiation from a straight segment 

A 

above a perfectly conducting plane driven by the propagating current pulse I(t - £j *r|/Vj). In the fre - 

jlcwE- *7 

quency domain this current is l(f)e 1 1 where rjj = c/vj. The solution is similar to that from a 
linear travelling wave antenna and can be obtained using a fraunhofer approximation (e.g., Le Vine 
and Meneghini, 1978a, b). One obtains 


E 7h,vi,z ( ^ * nAm/c) I(»0 ^ [ei kr+ - eJ kr “] 


(4) 


where 


T± = 


*V = 


fu =* ■ 


R L -A 
- ± - ( 1 , + t- ?R) + r 0 

A A AAA 

(g - z) -(£ • z)(z • PR) 2 

*?-(?• PR) 

(? • 7R) £ • 7R) - (i * z) (z • 7R) 2 


v - d • 7R) 


(5a) 


(5b) 


(5c) 


In this expression L is the length of the segment and t q is a constant needed to keep track of the 

phase of the current pulse as it propagates from element to element. The two terms in the solution 

corresponding to r± mathematically appears as if they come from the ends of the filament (e.g., see 

Le Vine and Meneghini, 1978a, b): The terms ±— ($ • VR) are the linear terms in a binomial expan- 

2c 

sion of the time for propagation from each of the ends of the filament to the observer and Lrj/2c 
accounts for the time (L/v) required for the pulse to propagate from one end of the filament to the 
other. The total electric field at the observation point is the sum over all elements. Using Equa- 
tions 2 and 4 it can be written in the form : 


W&j-M 





( 6 ) 


The voltage output of the radio receiver can now be computed by substituting Equation 6 into 


Equation 1. One obtains: 


(7) 


'JtCv. *i 


't 






V„(D ■ 


a(^ 0 ) K^q) f 

- ff ^_oo an i K i 1 1 1 


. , W.'tt',- _ c i2^^ + {)r,- ]H({)c -j:»ft d£ 


where it has been assumed that over the passband of the filter H(J), the factors a(v). I(i/) and the 
Ah jV (i>) are essentially constant. This is a good approximation for receivers such as those used by 
Le Vine (Le Vine et al , 1 976) and Le Vine and Krider ( 1 977) to measure time delay which had an 
effective bandwidth of about 300 kHz and made measurements at 3 MHz or more. (At 3 MHz the 
bandwidth is 10% and at 300 MHz the bandwidth is only 0.1%.) Recognizing that the integral in 
Equation 7 is a Fourier transform, one obtains: 


v Q (t) - 


V/e 




aO'o) H^) n 

2* afli K i 1 ill 


*(h(t-T i +)e^ ,r * ; o r i - h(t-r r )e^ 27riJ o T rj 


(8) 


where h(t) is the temporal, complex analytic, response of the filter (e.g., Bom and Wolf, 1959). 
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IV. RESULTS 


Calculations have been made of the receiver output, V 0 (t), for several different frequencies v Q 
using Equation 8 and letting the observer be on the x-axis 100 km from the base of the channel. 
The calculations were performed by first constructing a return stroke channel using linear segments 
placed end-to-end as described above. The change in cartesian coordinates of the end points were 
selected from a normal distribution with standard deviations o x = Oy = 3.8 m and o z - 6.0 m 
and means m x = m y = 0.05 m and m z = 8.0 m. The means of the x and y components were 
changed at randomly selected intervals using numbers from another normal distribution with zero 
mean and a x = o y = 0.58 m. These numbers were chosen to give agreement with statistics on the 
tortuosity of lightning channels obtained from photographs of lightning. The compound exponen- 
tial model I(t) = - e"^ 1 ) + I j(e^ yt — e - ^ 1 ) was chosen for the current. The current pulse 

is based on a form proposed by Uman ( 1 969) and others (e.g., Uman and McLain, 1969). The first 
two exponentials represent the main pulse in a form proposed by Bruce and Golde (1941) and the 
second pair of exponentials represent ..termediate current (Uman, 1969). The parameter values 
used were « = 4 x 10 4 sec' 1 , /? = 8 x 10 5 sec" 1 . y = 10 3 sec" 1 , 5 = 2 x 10 4 sec" 1 . I Q = 30 kA and 
1 1 = 2.5 kA which were based on the values suggested by Dennis and Pierce (1964) but modified 
slightly to reflect recent evidence that the rise time of the current pulse should be more rapid (e.g.. 
Clifford, Krider, and Uman, 1979; Weidman and Krider, 1980). The velocity of propagation of the 
current pulse up the channel was chosen to be 10 8 m/s which is representative of return strokes 
(Uman, 1969) and for convenience it was constant on all segments. The system filter was chosen to 
be h(t) = ( 1/T) sin 2 (rir/2T) 0 < t < T and zero otherwise. The form for h(t) was chosen largely 
for convenience of computation. It’s passband in the frequency domain is approximately 1/T wide 
which was chosen to be about 330 kHz (T ■= 3 x 10" 6 seconds) and the scale factor 1/T multiplying 
the sinusoid is used so that the filter has approximately unit gain over its passband. Finally, the 
factor a(v) which includes the detector gain and the effect of the antenna was chosen arbitrarily to 
be unity. 
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Figures 9 and 10 are examples of V Q (t) at 3 MHz (Figure 9) and at 130 MHz (Figure 10) from 
the channel shown in Figure 8. The lower trace in each figure is V Q (t) for a channel over a perfectly 
conducting ground (a * *>) and the upper example was obtained using the dielectric constant of soil 
given by Hoekstra and Delaney ( 1974; Figure 13). The time axes on these figures have been adjus- 
ted so that zero corresponds to the time when radiation from the beginning of the return stroke 
would first reach the observer. Notice that in the case of a perfectly conducting ground plane, the 
radiation begins abruptly at very nearly the time radiation is expected from the beginning of the re- 
turn stroke. A small delay is present because the system is causal and is essentially determined by the 
impulse response, h(t), of the filter. In contrast, when realistic values of e r are used for the ground 
plane (top curves in Figures 9 and 10) the receiver output shows a gradual it rease, building from a 
small value at the beginning of the return stroke to a function of time which eventually looks like 
that from a return stroke over a perfectly conducting ground plane. (Notice that the amplitude of 
V Q (t) as plotted in Figures 9 and 10 has been normalized by I(i>) the spectrum of the current pulse 
at the frequency, v Q . The normalization factors are given in the figure captions. Using the normali- 
zation factors one finds that the signal level at 3 MHz is about 1 .9 x 10' greater than at 130 MHz 
because I(i>) is very much smaller at 130 MHz than at 3 MHz. The same normalization has been 
used, independent of e r Hence from Figure 9 one sees that radiation at 3 MHz from a return stroke 
above a perfectly conducting ground earth is about 5 times as great as from the same channel over a 
more realistic soil.) 


A convenient interpretation of the results shown in Figures 9 and 10 can be obtained by notic- 
ing that Tj+ = r~n in Equation 8. This is true because the top of one segment coincides with the 
bottom of the adjoining segment. Using this result, Equation 8 can be written 


v 0 (t) - 




2ir 




( a i “ a i+i) 


all i 


R i 


, , j2jrp rt r 

h(t-fj+)e °i 


(9a) 


where 


(9b) 


a i = A h/ , 'o> e h i + \< v oK 1 
In this form, it is clear that the radiation detected by the receiver appears to come from the junc- 
tions between segments (Le Vine and Meneghini, 1978a, b) and that the receiver response to the 
radiation from such a junction is its impulse response, h(t), weighted by the factor (aj — a^jJ/Rj. 
This amplitude depends on the relative orientation of the segments, being large when the angle be- 
tween segments is large and small for segments which are nearly parallel (a^ = a^j). All junctions 
do not radiate simultaneously because the current pulse is propagating up the channel with finite 
velocity. Consequently, the sum in Equation 8 is only over those junctions which the current pulse 
has reached. Furthermore, h(t) has a finite width (i.e., is nonzero only in the interval t — T < t < 
t + T) and this further limits the junctions contributing to V Q (t) to those which have been illumi- 
nated within ±T seconds of t. Hence, one can imagine a window 2vT meters long moving up the 
channel at the velocity v of the propagation of the current pulse. The signal V Q (t) out of the re- 
ceiver comes from junctions which lie within this window. At the instant radiation begins (i.e., 
current starts up the channel) the receiver detects radiation from a single junction (the base of the 
channel) and the receiver output has the shape h(t). the receiver impulse response. As the current 
pulse advances up the channel radiation comes from an increasing number of junctions and the 
receiver output increases correspondingly. Eventually, a steady state is reached in which a new 
junction is added to the filter window t-T<t<t + Tby the advancing current pulse as old ones 
leave. The amplitude of V Q (t) in this “steady state” is determined by the values of a^ - aj + j. 
Finally, when the current pulse reaches the top of the channel, the amplitude decreases as radiation 
from fewer and fewer junctions is passed by the filter. This occurs at about 65 ps in Figures 9 and 
10 . 


When the surface is a perfectly conducting boundary Ajj = \ = 1 and only the relative ori- 
entation of adjoining segments affect the receiver output. In fact, a careful examination of Figures 
7-9 shows that the peaks in V Q (t) when a = correspond to elements that join at large angles. 
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When the boundary i$ not a perfect conductor, the interpretation of the receiver response as described 
above iz 5 ' ; ’> appropriate, but now the amplitude of the radiation is modified by the A^ v . For 
frequencies above a few MHz and observers far from the channel (i.e., R > z), one finds that A h a 
Ay - ^ 1 + f) where f = z'y'e^/R. Consequently, one expects V 0 (t) to have the same general struc- 
ture as in the perfectly conducting case but with its amplitude decreased by A v (t'). AyO) is a func- 
tion of the height of the radiating junction above the surface, and Figure 1 1 is a plot of Ay(p) as a 
function of height for several different frequencies. These curves were obtained using the values 
given by Hoekstra and Delaney (1974; Figure 13) for the dielectric constant of soil. One sees from 
Figure 1 1 that radiation from junctions close to the surface is severely attenuated and that this at- 
tenuation decreases with increasing height above the surface. At frequencies less than a few 1 00 kHz 
the attenuation is significant only very near the surface, whereas at frequencies of I MHz or more 
the attenuation is significant even at the top of a lightning channel. Because the current in return 
strokes propagates up the channel, height can be associated with time, the time it takes the current 
to reach a particular altitude. This time has been indicated on Figure 1 1 using v = 10 8 m/s for the 
velocity of propagation of the current pulse. For example, consider an observer at R = 50 km from 
the channel. Then the current pulse reaches 5 km (x = 0.1) in 50 ais(T = 1 .0). 

The effect on V Q (t) of this change in Ay(i») with time is clearly evident in the upper examples 
in Figures 9 and 10. In each case V Q (t) exhibits a gradual increase with time which is determined 
by the increase in Ay/p) that takes place as the current pulse propagates up the channel. Compari- 
son with the radiation when the channel is over a perfectly conducting ground plane (bottom of 
each figure) shows that V Q (t) has the same structure (i.e., same peaks and valleys) in both cases, as 
is especially evident at the ends of the record (e.g., near t = 50 ps), but with an amplitude which 
initially increases with time. This is so because the structure is determined by the shape of the 
channel which is the same in both cases and A(v) is changing from a constant value of unity for the 
perfectly conducting ground plane to a value which increases with height in the case of a realistic soil. 
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Another set of examples is shown in Figure 1 2 for a different lightning channel. The channel 
in this case is shown on the left and at the top right are V Q (t) for receivers tuned to 3 MHz and 130 
MHz and at the bottom are the receiver responses at these frequencies in the case of a perfectly con- 
ducting ground plane. The attenuation of the radiation emanating from the lower portions of the 
channel when the conductivity of the ground is taken into account is clearly evident in these exam- 
ples. 


Finally, Figure 13 shows the average receiver response obtained from three different channels 
all generated from the same statistical distribution. The attenuation of radiation coming from near 
the base of the channel is clearly evident when comparison is made with radiation from the same 
channel over a perfectly conducting ground plane (bottom). 
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V. CONCLUSION 


This paper has attempted to show that the time delay observed between arrival of fast electric 
field changes and RF radiation from first return strokes can be explained in terms of conventional 
models for return strokes by including the effect of the conductivity of the ground. This was done 
using a piecewise linear adaptation of the transmission line model for return strokes and using a 
solution patterned after that by Norton (1937) to account for the effect of the conductivity of the 
ground on the electric field radiated to a distant observer. The radiated electric field was then used 
as the input to a model AM radio receiver to calculate the RF radiation at different frequencies 
from the return stroke. It was shown that the radiation from different portions of the return stroke 
is attenuated compared to the radiation from an identical channel over a perfectly conducting 
ground plane. This attenuation is height dependent, being larger for portions of the return stroke 
near the base of the channel and decreasing for portions of the channel further from the surface. 
When this attenuation is combined with the fact that the return stroke current propagates up the 
channel at finite velocity, the result is a time dependent attenuation of the receiver output voltage. 
It was shown that the receiver output voltage is initially small because of severe attenuation of the 
RF radiation from elements of the channel near the ground and then the signal increases as the radi- 
ating elements are located further and further from the ground. This means that the receiver will 
not detect radiation from the return stroke until the current has propagated sufficiently high above 
the ground that radiation from this portion of the channel exceeds the detection threshold. If the 
conductivity of the earth were constant as a function of frequency, a similar delay would be expec- 
ted in the output from a conventional fast field change s> item. However, the conductivity of the 
earth is very strongly a function of frequency below a few Mi> and the imaginary part of the rela- 
tive dielectric constant of slightly moist soil is on the order of 10 5 at 1 kHz and 10 3 at 100 kHz. 
Consequently, at the frequencies which determine the gross shape of fast electric field changes it 
seems reasonable to assume that radiation takes place over a perfectly conducting ground plane. 
Thus, the fast field changes arrive at the observer with its shape largely unaffected by the ground 
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plane (of course the rise time and other high frequency detail will be affected) whereas the RF radi- 
ation, especially at VHF and above, will have been attenuated during the initial portions of the 
return stroke. The result is an apparent difference in the time of arrival at the observer between the 
beginning of the return stroke as deduced from the fast field change and the detection of the RF 
radiation. 
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Figure 1. Schematic of experiment to measure time delay. 
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Figure 2. Fast electric field change (bottom) and vertically 
polarized radiation at 295 MHz (top) from lightning in 
Florida at the Kennedy Space Center. The data was recorded 
during TRIP-76 on July 16. 1976. 



Figure 3. Fast electric field change (bottom) and horizontally 
polarized radiation at 295 MHz from lightning (top) recorded 
at the Kennedy Space Center, Florida. The data were recorded 
during TRIP-76 on July 10, 1976. 
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Figure 4. Fast electric field change (bottom) and vertically 
polarized radiation at 3 MHz from lightning in Virginia The 
data were recorded on July 29, 1979 at the Wallops Flight 
Facility, Virginia. 



Figure 5. Fast electric field change (bottom) and vertically 
polarized radiation at 139 MHz from lightning in Virginia. The 
data were r fc . ,rded on July 29, 1979 at the Wallops Flight 
Facility, Virginia. 
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Figure 6. Fast electric field change (bottom) and vertically 
polarized radiation at 3 MHz (top) from lightning in New 
Mexico. The data were recorded near Langmuir Laboratory 
on July 21.1979 during TRIP-79. 



Figure 7. Fast electric field change (bottom) and vertically 
polarized radiation at 1 70 MHz (top) from lightning in New 
Mexico. This example was recorded during TRIP-79 near 
Langmuir Laboratory on July 21, 19" T 9. 
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Figure 8. A lightning channel constructed by placing straight segments end-to-end in a random 
walk. There are approximately 450 segments in this example. 
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Figure 10. The receiver output at 130 MHz for radiation from the channel shown in Figure 7. The 
oh*»Tv a r is 100 km from the channel on the ground. Multiply the receiver output by 3.42 x !0“ 8 

to get V 0 (t) in volts. 
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Figure II. The correction factor, A v (v) which accounts for the conductivity of the ground 
plane on radiation from a dipole, h is the height of the dipole and R is the distance to the 

observer. 



















APPENDIX A 

MODEL FOR AM RADIO RECEIVER RESPONSE 


It is the purpose of this appendix to derive an expression for the voltage response of an RF 
receiving system such as used in the experiments to measure the time delay. The RF receiving sys- 
tem used in these experiments consisted of an antenna, usually vertically polarized, connected to a 
standard AM radio receiver and followed with a post detection filter. The post detection filter was 
employed to reduce all measurements to a common 300 kHz bandwidth. 

The antennas used were vertical dipoles or vertically polarized disk-cone antennas. These are 
isotropic antennas which, over the bandwidth of the measurements, can be assumed to deliver a 
voltage proportional to the incident vertically polarized electric field. Except for a phase which is 
ultimately lost in the detector, this proportionality constant isi/jrG/A e /k where k = a\/p7 = 

2 ir/X, G is the gain of the antenna in the horizontal plane (parallel to the ground), and A e is the 
effective receiving area of the antenna (Kraus, 1966). 

The AM receiver is a device for detecting and amplifying the envelope of an amplitude modu- 
lated sinusoid (carrier) at a particular frequency. This is normally done by translating the input 
signal to an intermediate frequency where the actual processing is done. However, the frequency 
translation is done for engineering purposes to make the detection more efficient, and it is not 
necessary to do this in order to model the receiver output. The ideal device is a perfect envelope 
detector in series with a filter which represents the equivalent bandwidth of the system and an ideal 
amplifier which represents the system gain. 

For systems whose bandwidth 6 is small compared to the nominal frequency, v Q , of the mea- 
surement, these operations can be written explicitly in terms of the spectrum (Fourier transform) 
of the incident radiation. To do so it is convenient to write the vertical component of incident 
radiation, E z (t), in the form: 


A-l 
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E z (t) = 2Re r E z O)e~i 2ff,;t dv 


(Al) 


where Re means “real part of’, where E z O) is the Fourier transform of E z (t), and where the inte- 
gral and factor 2 are called the complex analytic representation of E z (t) (e.g., Bom and Wolf, 

1959). Using this notation, the signal V(t) out of the antenna and filter and into the envelope de- 
tector is: 

V(t) = 2Re f aaOHME-Me-J 2 **' 1 dv (A2) 

J 0 

where HO) is the Fourier transform of the filter h(t) and aO) “ (A/k\/jrG/A e is the combined 
effect of the antenna and an amplifier with gain, A. Although in practice the post detection filter 
is applied to the video output of the receiver, mathematically the effect of the post detection filter 
can also be included in HO)- This will be done here and it will be assumed that the equivalent filter 
HO) has a passband centered about frequency v Q which is very narrow compared to v Q . That is: 

HO - v 0 ) 


HO) = 


- v Q I < 6 


(A3) 


\v - t> 0 l> 6 


and 6,0 0 < 1. Using this assumption and making the change of variables £ = v - v Q , Equation A2 
becomes 

V(t) = Re 


where 


- e -j2^ 0 t f 2 aO 0 + £) E z 0 0 + £) H(£)e“j 2ir ^ d£ 
= I e z 0 0 .t) I cos(27rv 0 t + 0) 

yt oo 

e z 0 0 ,t) = / a0 0 + 1) E z 0 0 + £) H (£) e-i 27 * d£ 

AM 


(A4) 


(A5) 


and where the assumption that HO) is nonzero only in a narrow band about v Q (Equation A3) has 
been used to formally extend the lower limit of the integration in Equation A4 to From Equa- 
tion A4 it is clear that V(t) has the form of a sinusoid, cos(2irv 0 t + 0), at frequency i> 0 , which is 
amplitude modulated. The output of the detector is the envelope, I e z O D ,t) I, of this carrier. Thus, 


the output. V Q (t). can be written: 


ORIGINAL PAGE 53 
OF POOR QUALITY 


V 0 (t) 



2a(«/ 0 + *) E z (u 0 + *) H(t)e-i 27r t* d* 


(A6) 


which is the form used in the text to compute the receiver response to an input signal, E z (t). 
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APPENDIX B 

RADIATION FROM A CURRENT FILAMENT ABOVE A CONDUCTING EARTH 


The building block used in this paper to construe* a lightnii.* channel is a straight filament on 
which it is assumed that current propa* ues at constant velocity. The purpose of this appendix is to 
present a solution for the radiation from such a filament which is valid when the filament is over a 
finitely conducting ground plane. A solution is desired for the vertically polarized components of 
electric field (i.e., z-component) at an observation point which is on the surface (z = 0) and far from 
the lightning channel compared to the height of the filament above the ground plane. 

The desired solution can be constructed from the solution for radiation from a dipole (i.e., a 
point source) with polarization F, above a conducting ground plane as illustrated in Figure B 1 . 
Radiation from the dipole is a classic problem in electromagnetic field theory originally treated by 
Sommerfeld ( 1 926) and Norton (1937) and in other forms by Wait (1961) and Banos (1966). The 
derivation presented here is after that by Norton (1937). 

Since radiation from a dipole is desired, it is convenient to solve for the fields in terms of a 
Hertz potential (e.g., Stratton, 1941). One obtains: 

B = -ju>A/ 0 e 0 e r ^ x n (Bl) 

E = k 0 2 e r ff + V(Vxir) (B2) 

where the Hertz potential it is a solution to the wave equation 

V 2 it + k 0 2 e r 7r = -P/e 0 e r (B3) 

and k 0 = gVm 0 c o and e r = € / f o + -i a / w€ o = e r + J e r "• 1° the region above the conducting plane 
(z > 0) it will be assumed that e T = 1 . 

WhenP = P^ (that is, when the dipole is vertical) only one component v = ir^z is needed 
and one can show (e.g., Banos, 1966; Norton, 1937) that 

r v z = [gfT/n-gfT/Tj') + V(T/r I ')] (B4) 

€ o 


B-l 


where 


V( r/ r j') = r— f — J 0 (Xp)e’ (Z +Z)8 X d\ (B6) 

1 2tr ./q m + ® c r 

andp * v/(x.'- x)^ + (y'~ y) 2 , £ 2 = X 2 -k 0 2 , m 2 = X 2 -k Q 2 e r and where «( 7/7') = 
exp0k o R)/(4?rR) where R 88 17 -7' I is the Green’s function for the wave equation in the absence of 
boundaries (i.e., the free space Green’s function). Also 7, ' is the position the mirror image of the 
dipole below the ground plane and I 7 -7j'l * Rj is the distance between the observer and 
this image (Figure Bl): 7,' = 7' - 2(7'*$)$. Using Equation B 7 the following form is obtained 
for the vertical component of electric field from the vertical dipole in the region z > 0: 

E v z = [k Q 2 + 0 2 /3z 2 )] (B6) 

The solution when the dipole is horizontal requires two components of the Hertz potential. 

For example in the case of a horizontal dipole oriented along the x-axis (P - P x x) the solution 
can be expressed in terms of an x- and z -component of if : iF* 1 = z. However these 

components are related and the solution for E z in this case can be written in terms of the integral 
V( 7 /7»'), just as in ihe case of a vertical dipole. Following Norton (1937) one obtains: 

. p x & 

gh _ C os0 — - [V(r/r,') - g( r / r ') - g( r/ rj')] (B7) 

z e 0 a P dz 

where p is defined as before and p cos 0 88 x - x\ In this notation, 0 is the conventional cylindri- 
cal coordinate measured from the x-axis. In the case of a horizontal dipole along the y-axis, the 
only change is that cos 0 is replaced by sin 0. For the application to the problem of the time delay 
in radiation from lightning, it is desired to evaluate Equation B6 and B7 when kR is large. In this 
case the following asymptotic approximation for V(7/7') pertains (Norton, 1937): 

V(7/7') as 2*7/7,') (1 - */(!+$)] (B8) 


where 


B-2 


* 
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z + z ' 
R I 




(B9) 


Substituting the approximation into Equation B6-7 and letting z * 0 (observer on the surface) one 
obtains 

E; - -2w 2 M 0 P z [f/n+f)]g(T/r') (BIO) 

R 


E£ = -2w 2 *i 0 [$/(1+$)] 1 +j 


k 0 (z') 2 (I+$) 


(z7R)(P h * PR)g( 7/7') (Bll) 


where P^ = Pjj X x + P^y y. As a check on these results, notice that in the limit e r j°° one obtains: 

(B12) 


JkR 

E z jcoAt 0 (2I v ) = Eyz 


eJ kR 

E z - jwM 0 (2 l h • ^R) (z'/R) — = E£ (B13) 


where I = -j«P is the current equivalent of the time varying dipole and where E^ and E h"z are 
the standard results one obtains for radiation from a dipole above a perfectly conducting half-space 
(e.g.. Jones. 1964). 


Comparing B 1 0-1 1 with B 1 2— J 3 it is clear that one can write 

E; = ^E” (B14) 

Ej = A„E h ~ (B15) 

where 

* 

\ " 7 TTf (B16a) 

* = z'v^/R (B 1 6c> 

and z' is the height of the dipole above the surface. 


B-3 
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It is interesting to compare the relative contributions of vertical and horizontal dipoles to E z . 
Thus, from Equations B 1 2-1 6 and letting £ < 1 , one obtains 



(B17) 


Hence, if the dipole is several wavelengths above the surface, then for an observer far from the di- 
pole compared to its distance above the surface (R > z), the radiation will be much stronger from 
vertical dipoles than from horizontal dipoles of equal strength. In the application of the results to 
lightning return strokes one would expect the vertical components to be even more dominant be- 
cause the channels have a predominant vertical orientation. In fact, it is a common assumption in 
calculating radiation from return strokes to assume that the channel is straight and vertical (e.g., 
Lin, Uman and Standler, 1980). 


The analysis to this point has been for radiation from dipoles (i.e., infinitessimally small, point 
sources). The solution that is needed for the return stroke is for radiation from a filament (straight 

- A /\ 

line) along which current propagates at constant velocity: That is, I(y,T') = £ I(i>) exp(jk 0 Tf£ * 7') 

A 

wnere £ is a unit vector along the filament in the direction of propagation of the current and rj = 
c/v is the ratio of the speed of light in vacuum to the velocity of propagation, v, of the current pulse 

A. A 

along the filament. In the time domain this current pulse is I(t) = £ I(t - £• r /v). The solution 
for the dipole presentea above is in effect a Green’s function for this problem and a formal solution 
for the z-component of radiated electric field can be written in terms of this solution by summing 
along the filament. One obtains: 


B-4 
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where 


^(h.v)z 

gh( 7 / 7 ') 
gy( T/7*) 
lO'.r') 


/ 


filament 


gh, v <f/T') !(*>/) dr 


^ ^ gjkoR 
}UH 0 [t-z) ;ffR \ 

/v ~~ ^ e jk ° R 

ju/i 0 l(C* PR) - (£*z) (z* PR)] (z'/R) jirK A h 


(B 18 ) 

(B19a) 

(B19b) 

(B19c) 


In the limit Ay = = 1 the Green’s functions v (7/7') are just the Green’s functions for the 

z-component of electric fielc adiated from a current above a perfectly conducting ground plane. 
Denoting these Green’s functions by g£° v (7/7') and substituting Equation B19c explicity into 
Equation B 18. one obtains: 

E (h,v)z * J ^h.**h!v(T/T') I W e ,k °’' <! ' r d7' (B:0) 

filament 

For filaments which are sufficiently short that £ can be assumed to be constant, the A^ v can be 
taken out of the integral. After doing so, the remaining integral in Equation B20 is the conven- 
tional solution for the fields radiated to an observer on the ground from a linear travelling wave 
antenna above a perfectly conducting half-space. Using a fraunhofer approximation to evaluate the 
integral, which is valid when the observer is far from the filament compared to its length (e.g., 

Le Vine and Meneghini. 1978a, b), leads to Equations 4 and 5 in the text. 
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Figure B1 . The geometry for calculation of radiation from a dipole over a conducting 

ground plane (2 < 0). 


